B —> Dru Branching Ratios: 
Opportunity for Lattice QCD and Hadron Colliders 
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<D ' Abstract 

^sO \ In the Standard Model, scalar contributions to leptonic and semileptonic decays are helicity 
suppressed. The hypothesis of additional physical neutral/charged Higgses can enhance such 
scalar contributions and give detectable effects especially in B physics. For the charged Higgs, 
experimental information on both Br(B — ► Dtv) and Br(B — > tv) has already become available 
and in particular the B — > Dtv branching ratio measurements will be further improved in the 
coming years. Hadronic uncertainties of scalar contributions in semileptonic decays are already 
in much better shape than the ones plaguing the helicity suppressed leptonic decays B — > tv. 
Combining existing experimental information form the B factories, we explore which existing and 
future lattice estimates will be useful to directly address new physics effects from measurements 
, of Br(B Ui d iS — ► D u d s Tv), which can be performed also at hadron colliders. 
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As is often stressed, in the near future the LHC will represent the main avenue to establish 
00 ! the presence of new physics by directly detecting new particles at the TeV scale. On the other 
hand, virtual effects of these particles can affect low-energy observables, probed mainly by the 
flavor factories and soon by the LHCb. As has been proven by the B-factories, the energy reach 
of such indirect searches can often surpass direct detection strategies, making them worthy of 
' pursuit even at the opening of the new energy frontier. Among the possible new particles, 
the Higgs boson is the only one expected in the Standard Model (SM) picture. At the same 
time, we have to observe that the established SM parametrization of the Higgs sector is only 
a conservative example of a possible electro-weak symmetry breaking mechanism. The present 
information on the massive W and Z bosons from electro-weak precision tests only constrain 
the goldstone modes p] of the Higgs field while leaving space for an extended physical Higgs 
sector. Namely, additional neutral/charged Higgses appear in many models trying to solve the 
inconsistencies of the SM. 

Therefore theoretical and experimental study of scalar effects in observables, mediated at 
tree- level by neutral/charged boson^] is vitally important in future experimental programmes. 
In particular, effective density operators from charged scalar boson interactions have to be 
considered in the effective weak Hamiltonian, which for b — ► q(u, c) transitions, for example, 
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1 Loop induced flavor changing neutral current processes, for example [2] b — » sj, can be sensitive to additional 
Higgses but this information is diluted by contributions from other particles and final constraints are model- 
dependent [3]. 
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In the minimal flavor violating (MFV) extensions of the SM [3] by an additional Higgs doublet 
the additional new physics (NP) coupling can be written as 



C 



m b m,£ tan 2 /? 



NP 



m 



2 H+ 1 + e tan (3 ' 



(2) 



where tan (3 is the ratio of the two Higgs VEVs while eo parameterizes possible Peccei-Quinn 
symmetry breaking corrections and is typically of the order of 1% in the MFV minimal su- 
persymmetric SM (MSSM). Due to the suppression of quark and lepton Yukawa couplings in 
eq. ([2"]). B helicity-suppressed processes receive largest effects from the charged Higgs. In this 
respect, the B — > tu decay branching ratio [5], given by 
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has often been stressed as a good candidate and the recent B-factory results have given important 
constraints on CJ^ P . Unfortunately, the presently established experimental precision is only 
about 30% and unlikely to improve in the near future as the perspectives to measure B — > tu 
at the Tevatron or LHCb are highly compromised. Furthermore, the SM expectation estimate 
presently suffers from sizable parametrical uncertainties induced by \V ub \ and f P - This opens 
the door for alternative modes to be studied with the present experiments. 

While Higgs effects in K and D modes are small and difficult to disentangle at present 
theoretical precision [61 [7] , the situation is much better in the case of semileptonic B — > Dtv 
decays [SI [10]. The partial rate can be written in terms of w = vb • vd as 
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where t(w) = m 2 B + m 2 D — 2wmr)'mB and we have decomposed the rate into the vector and scalar 
Dalitz density contributions 
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where G(w) and A(w) encode our ignorance of the QCD dynamics. Even before analyzing the 
theoretical uncertainties of these modes let us note that the present constraints on C^p from 
[6] and B — > tu [Tl] decayal still allow for sizable new physics effects in eq. @ for the 



2 In details, the eotan/3 terms in eq. ([2} are set to be equal between B — > tu and K — > /iv, as it happens in 
MFV MSSM. 
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□ Allowed region 
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Figure 1: In green we plot as a function of w the allowed region for p$ p (w)/ ps(w) in eq. (J7|), 
using constraints from both B — > tv and K ^ pv decays |1 1|. [6]. A large deviation from the 
unity, the SM expectation, is still possible with respect to the SM. Note that ps{w) contributes 
50% to the Br(B -» Dru). 



case of B — > Dtv as represented in fig. [TJ where the allowed region of the helicity suppressed 
contribution of eq. ([!]) for B — > Dtv, namely 
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is shown. 

The main parametric uncertainties in eq. (j3J) are represented by the modulus of V c b and the 
hadronic form factors G(w) and A(w). Presently, the most accurate value of \V c b\ = 4.15(7)% 
comes from the fit to inclusive B — > X c iv decays which are insensitive to scalar contributions [9] . 
Because of charm states, information from heavy quark expansion for the form factors is a 
priory unsatisfactory, since corrections to the static limit m c ,nib — > oo, formally parametrized 
by £ = l/mfe(l — mb/rric) can be large and undetermined. More reliable information is expected 
from the lattice and indeed a number of studies have computed the normalization of the vector 
form factor G(w) at to = 1 to a precision of a few percent, while a recent study extended its 
determination to a region of w G [1, 1.2] [12] . These values must then however be extrapolated 
over the entire kinematically accessible decay phase space, which is larger in the case of B —> Deu 
(w £ [1, 1.59]) than for the tau mode (w £ [1, 1.43]). For such an extrapolation, HFAG adopts 
the parametrization of G(w) [13] H 

G(w) = G(l) x [1 - 8p 2 z{w) + (51p 2 - W)z{w) 2 - (252/> 2 - 84)z{wf] , (8) 



3 Using analyticity and crossing symmetry, a general parametrization for semileptonic decays has been proposed 
in ref. Q3] . However, for modes such as B — > D, the smallness of z, and the judicious use of heavy-quark symmetry 
in ref. [13| . allows for a especially tailored parametrization in terms of eq. ((8} 
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Figure 2: Comparison of B — ► Dez^ form factor determination from Belle [T7], Cleo [18] and 
Lattice QCD. The latter data points have been multiplied by the HFAG world average value of 
\Vdt\ from inclusive measurements. The HFAG average fit [16] to eq. (JSj) is also shown. 



with z(w) = (y/w + 1 — y/2) /{\J w + 1 + y/2) in terms of two parameters: the normalization 
G(l) and the slope p 2 . In addition, in the SM and as well as in its MFV extensions only 
G{w) will actually contribute to B — > Deu and one can use experimental information on the 
differential decay spectra in such an extrapolatiorQ. At present, the HFAG [16] experimental 
information consists of relatively old publications by Belle |17] and Cleo [18]. We can use this 
information however to asses the relative precision obtainable from combining lattice information 
with experimental inputs efficiently We compare in fig. [2j the Belle [17] and Cleo [18] data on 
\V cb G(w)\ and the HFAG fit to the data from eq. © (using |T^,|G(1) = (42.3 ± 4.5)1(T 3 , 
p 2 = 1.17 ±0.18 with correlation 0.93), together with the lattice data from ref. [12] and the 
fit from eq. © to the lattice results of G(w) p2] (yielding G(l) = 1.03(1), p 2 = 0.97(14)) 
both multiplied by the HFAG value of \Vd,\ mentioned above. The two sets are in agreement 
at present precision (10% on the normalization and 15% on the slope). Improvement however 
could come from several sources: Babar has already announced to improve the measurement of 
the differential decay rate to allow for extraction of p 2 to below 10% by reducing the statistics 
error of Belle by a factor of 4 [19jl. However, to be able to apply this precision to the integrated 
rates, one would need to precisely determine either G(w) on the lattice while using inclusive 
determination of |V^,| or consider ratios, where the overall normalization factors of |V^,G(1)| 
cancel. 

On the hand, the uncertainties coming from A(u>), which regulates the helicity suppressed 
terms, are already much smaller, especially than those plaguing the dimensional variable Jb 



4 For completeness, the mechanism introduced in ref. [15] to enhance electronic modes in K — > ev T and 
B — > ev T by orders of magnitude gives negligible effects less than 0.1% for the partial rate of B — > D£v, once the 
K — > ev T bound [6] is taken into account. 

5 At this level of precision, non-helicity suppressed NP contributions to the b — > ceu transition could be 
constrained for the first time (for example R-parity violating MSSM [20]). 
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entering B — > Iv decays. In other words, the current (quenched) lattice estimate of A(w) for w 
in the range 1 — 1.2 is at about 2% precision, consistent with a constant value of A(w) = 0.46(1). 
Mainly, such an achievement on the lattice was possible by introducing double ratios of lattice 
correlators [21] and 6 boundary conditions [22]. Moreover, this precision can further be improved 
by studies involving unquenched simulations and lighter sea quark masses. In particular, a 
measurement of B s — > D s £v will opt for lattice data on B s — ► D s form factors including scalar 
contributions. These however no longer require chiral extrapolations for the valence quarks, 
eliminating important sources of systematics. Finally, since A(w) only contributes significantly 
to the decays involving taus, the extrapolation from the region presently probed by lattice 
simulations to the complete kinematically accessible region is not large as is the case for the 
G(w) form factor in B — > Dev transitions. 

We finally combine these lessons and try to project the present sensitivity of B — > Dlv decays 
to scalar contributions into the near future. We start with the ratio Br(B — > Dtv) / Br{B — > 
Dev) [9] [TO] which, as stressed above, even in the presence of NP scalar contributions only 
depends on two hadronic parameters, the precision of which can furthermore be improved in 
the near future: p 2 and A(w). By integrating eq. (j3]) with the use of eq. ([8]), the fitted lattice 
results for the form factor, A(w) and the HFAG value of p 2 as determined from the B — ► Dev 
spectrum, we average over the B& jU — ► Da, u modes to obtain 

bI{bZd7J) = ( °- 28 ± °-° 2) X [1 + ^WMChp) + 0.88(2)|C^| 2 ] . (9) 

We see that the SM prediction uncertainty is already below 8% and is expected to be improved 
soon with the new Babar data on p 2 . Interestingly, Babar has already published a value [23] 
for the above ratio with uncertainties of 30%, making it possible to compare with the B — > tv 
measurement and its bound on C r NP in fig. [3j Even more importantly, unlike B — > tv, this 
measurement can be improved at hadron colliders together with B s — ► D s tv. Therefore we plot 
the present exclusion region in the tan (3 - m H + plane in fig. H] together with the percentage 
deviation from the SM prediction for Br{B — > Drv)/ Br(B — > Dev) in the presently allowed 
region. 

An even more prospective observable however, may be represented by the ratio of partial 
B — ► Dr(e)v decay widths integrated over the same kinematical w region. Since in the case of 
B — > Dtv the kinematically available region is much smaller than for the B — > Dev one can just 
consider the full Br(B — > Dtv) |24j . while imposing a kinematical cut of w < 1.43 in the light lep- 
ton case. In this way one avoids the large extrapolation away from the lattice data points and fur- 
ther reduces the uncertainty due to the p 2 parameter. Presently such a ratio can be estimated at 
Br(B -» Dtv) / Br(B -» Dev)\ w<lA3 = (0.56 ± 0.02) x [l + 1.38(3)i?e(C^ P ) + 0.88(2) |C^ P | 2 ] 
with an error on the SM value of only 4%, while the relative new physics contributions are not 
affected by the cut at all, since they only appear in the tau mode. Once the experimental preci- 
sion for this observable would approach the above theoretical errors, one could further restrict 
the kinematical region considered closer to the one accessible to the lattice studies or finally 
consider binned or differential rates. 

In existing literature, the differential rates [9] [10] are often stressed as being highly sensitive 
to scalar contributions in B — > D transitions compared to the integrated rate. However such 
measurements will only become available with the advent of the Super Flavor Factories, where 
both the B — > Dev and B — ► Dtv spectra will be available at a few percent level in several 
w bins. Then, measuring the ratio of B — > Dtv and B — > Dev differential distributions [9] 
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Exp. vs Theory at 95% C.L 




Figure 3: The ratio Br(B — > Dtu) / Br{B — > Dez/) is shown together with the Br{B — ► rz/) 
as function of CJ^ P , eq. ([2]). Both curves have been normalized to their SM central values. 
Error bands on the curves represent the theoretical uncertainties at 63% and at 95% C.L.. The 
horizontal bands represent the corresponding experimental values 123]. 
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Figure 4: Exclusion region in the — tan (5 plane due to present determination of B — > tv 
(in blue) and Br{B — ► Dtv) / Br(B — > Deu) (in gray). Note that the small allowed band in 
the middle is excluded by K — > \xv determination [6] (not shown). Red dashed lines represent 
percentage deviation from the SM prediction of R = Br(B — > Dtv) / Br{B — > Deu) in the 
presently allowed region. 
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Figure 5: The quantity pg P (w) from eq. is shown for three values of w as a function of 
C7VP' ec l- ©• The va l ues of io = 1.1, 1.2, 1.3 are chosen to coincide with the presently available 
lattice data [12] . Experimentally, pg P (w) can be accessed at a Super Flavour Factory via the 
measurement of dT(B — ► Dtv) / dY{B — > D£u), eq. (j4|) in those u>-bins. 



integrated over given u;-bins gives direct access to Pg P (w) which can be compared with the 
lattice estimates of ps(w) in the same bins to obtain bounds on C T NP by reducing ambiguities 
due to G(w) estimates and w parameterizations. We project the potentialities of measuring 
Pg P (w) in eq. ((7|) with respect to Br(B — > rv) in fig. EJ 

In the meantime, the ratio of (partialy) integrated rates Br{B — > Dtv) / Br{B — > Dev) seems 
to represent the best strategy for indirectly probing charged Higgs contributions to low energy 
observables at the Tevatron and LHCb. Even if Br(B — ► Dtv) / Br(B — > Dev)\ w< iA3 can not 
be measured directly, precise data on V c b and B — > Dev decay spectra from the B factories can 
be used to obtain comparable precision directly on the B — > Dtv branching ratio. Moreover, 
since the bounds from B — > tv are affected by larger theoretical uncertainties, the B — > Dtv 
modes allow for an important crosscheck. Let's mention that at 95% with the present central 
value and with a smaller experimental error of 20%, the exclusion region from B — ► Dtv is 
already competitive to the one from B — > tv, while at 5% error the SM and the MFV MSSM 
would actually be excluded. Thus such a precise measurement of Br(B — > Dtv) together with 
further lattice studies of G(w) away from w = 1 and A(w) would be highly welcome since both 
the central values as well as an accurate estimation of their errors are essential to obtain valid 
bounds on new physics. 
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